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The thermal conductivity ,k, of Bi2Sr2CaCu208 was studied as a function of magnetic field. Above 
5 K, after an initial decrease, k(H) presents a kink followed by a plateau, as recently reported 
by Krishana et al.. By contrast, below IK, the thermal conductivity was found to increase with 
increasing field. This behavior is indicative of a finite density of states and is not compatible with the 
existence of a field-induced fully gapped d x 2_ y 2 +id xy state which was recently proposed to describe 
the plateau regime. Our low-temperature results are in agreement with recent works predicting a 
field-induced enhancement of thermal conductivity by Doppler shift of quasi-particle spectrum. 
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Until very recently, thermal conductivity in the vortex 
state of high-T c cuprates was analyzed through a picture 
which was originally developped in the context of con- 
ventional superconductors. According to this scheme, 
vortices constitute new scattering centers for heat car- 
riers and their introduction leads to a decrease in heat 
conductivity. Even in the absence of a satisfactory ac- 
count of vortex-quasi-particle interaction in unconven- 
tional superconductors, this picture provided a qualita- 
tive explanation of the observed field-induced decrease 
in k(H) in the investigated region of field-temperature 
plane . However, recent results reported by the Princeton 
group on the thermal conductivity in the mixed state of 
Bi 2 Sr 2 CaCu20 8 (Bi2212) § and La 2 _ x Sr x Cu0 4 @ have 
demonstrated the need for a more vigorous exploration 
of transport in the vortex state of the cuprates. Notably, 
their study of k{H) in Bi2212 |jj indicated that for tem- 
peratures below 20K, above a threshold field Hfc(T), k 
becomes insensitive to the magnitude of magnetic field 
applied along the c-axis. The authors proposed that 
Hk — which presents a rough T 2 dependence — is the sig- 
nature of a phase transition and the high-field plateau 
regime represents a new quasi-particle- free superconduct- 
ing state. In their suggested scheme, the low-energy exci- 
tations of the parent d x 2_ y 2 order parameter is supressed 
at Hk(T) due to the opening of a large d x 2_ y 2 + id xy gap 
H in the whole Fermi surface. 

In this letter we present a study of the field-dependence 
of thermal conductivity at various temperatures down to 
O.IK. Our results reveal an increasing k(H) at subkelvin 
temperatures. We will argue that this is a strong indica- 
tion for presence of field-induced delocalized excitations 
and as such rules out the hypothetical fully-gapped state 
invoked jl],[| to explain the plateau observed at higher 
temperatures. Moreover, our results highlight the rele- 
vance of field-induced Doppler shift of energy spectrum 
which — as recently pointed out [Q — has been hitherto ne- 
glected in analyses of thermal transport in the cuprates. 

Fig. 1 presents the thermal conductivity of two 



optimally-dopped Bi2212 crystals at temperatures just 
above 5K (T c is 89K for sample 1 and 88K for sample 
2). Our results in this temperature range confirm the 
behavior originally reported by Krishana et al. [ IJ : after 
an initial drop, k(H) presents a kink and then becomes 
quasi-constant. These authors attributed the plateau of 
k(H) to a sudden disappearance of quasi-particle (QP) 
transport, presumably consequent to the opening of a gap 
induced by the magnetic field. However, as we already 
argued in a comment presenting the data on sample 1 [g| , 
the existence of a strong hysteretic behavior weakens se- 
riously this interpretation. Indeed, if the high-field state 
was a superconducting state free of quasi-particles , one 
can not see why it should show such an extreme sensitiv- 
ity to the details of the sample's thermo-magnetic history. 
In their reply ||, Krishana et al. confirmed the exis- 
tence of a hysteretic behavior with a sample-dependent 
intensity. They suggested that the disorder in the crys- 
tallographic structure could induce textures and defects 
in the vortex lattice forbiding the phase transition oc- 
curing at Hk{T) to go to completion everywhere in the 
crystal. In this way, d x 2_ y 2 domains with low-energy 
excitations would survive in the gapped d x 2_ y 2 + id xy 
superconducting state. They further argued that in the 
sweep-down trace of k(H), the flux distribution would be 
more uniform, with a lower density of d x 2_ y 2 domains, 
and consequently, a lower density of electronic heat car- 
riers. However, this proposed explanation of hysteresis is 
unsatisfactory. It supposes that sweeping the magnetic 
field up to lO-fffc(T) leaves the density of the hypotheti- 
cal domains unchanged — which is necessary to keep a so 
flat plateau. Moreover, it fails to explain the sign of the 
hysteresis in sample 2. 

To follow the temperature dependence of Hk{T) we 
studied the thermal conductivity of sample 2 at lower 
temperatures (T j IK). For these measurements, instead 
of sweeping the magnetic field at a constant temperature, 
we used an alternative procedure: a constant magnetic 
field was applied above T c and then thermal conductivity 
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was measured as a function of temperature. This was to 
avoid the well-known problems related to inhomogeneous 
penetration of vortices in a zero-field-cooled sample. In- 
deed, by studying the magnetization of the sample at T = 
5K, we found the field of complete penetration to be 3500 
Oe. For this reason, it is not possible to resolve a possi- 
ble kink in k(H) for fields smaller than 3500 Oe which is 
the expected magnitude of Hk below 4 K. On the other 
hand, by cooling the sample below T c in a finite magnetic 
field, one would expect to attain a homogeneous field in 
the sample even at small fields. 




a field-induced increase in thermal conductivity is to in- 
voke an enhancement of the available excitations of the 
superfluid condensate. There are two possible origins for 
unpaired electrons in the mixed state. First, there are lo- 
calized excitations of the normal core. The second source 
is provided by the Doppler shift of QP energy spectrum 
due to the superfluid flow around the vortices. 




FIG. 2. Temperature dependence of the thermal conduc- 
tivity k(T)/T with a magnetic field applied above T c . Note 
the crossing of the 0.6 kOe and 0.0 kOe curves. 



FIG. 1. Field dependence of the thermal conductivity 
k(H)/k(H — 0) for two different samples. Note the initial 
decrease, the kink, the "plateau" and the large hysteresis. 

Fig. 2 shows the temperature dependence of thermal 
conductivity k(T)/T for various magnitudes of magnetic 
field. We found that even at the lowest temperature in- 
vestigated (T ~ 0.1K) k(T) does not present a cubic 
term characteristic of phonon conductivity with constant 
mean-free-path. Thus, the separation of the electronic 
and lattice components of zero-field thermal conductiv- 
ity is not straightforward. Here we will concentrate on 
the effect of the magnetic field which can be analyzed 
indepedently. 

As seen in the figure, at this temperture range, k is 
mostly enhanced when the sample is cooled down below 
T c in a finite magnetic field. Fig [| shows the field pro- 
file of thermal conductivity extracted from n(T) data at 
various fields together with n{H) at 5.5 K already shown 
in Fig. 1 and obtained by sweeping upwards the mag- 
netic field. The figure shows that at low temperatures, 
instead of presenting a plateau, k(H) increases with the 
magnetic field. Furthermore, this field-induced enhance- 
ment of thermal conductivity becomes more pronounced 
with decreasing temperatures. This behavior has strong 
implications for the debate on the origin of the plateau 
feature observed at higher temperatures. 

We begin by noting that lattice heat transport can only 
be reduced by magnetic field. The only way to explain 



Vortex cores can be excluded as a plausible source of 
heat transport in our context of investigation. Let us re- 
call that in a s-wave superconductor, the energy of the 
first bound state inside the vortex cores is A^/2Ef [Q, 
which for cuprates would yield a fraction of meV. In a 
d x 2_ y 2 type superconductor, according to analytical || 
and numerical (^] works, one do not expect bound states 
inside the vortex cores at all. Nevertheless, tunneling 
spectroscopy measurements have shown the existence of 
bound states inside the vortex cores of YBCO flipf , but 
the energy of the first one is 5.5 meV (60K). In Bi2212, 
the very recent detection of vortices by tunneling spec- 
troscopy jllj revealed a pseudo-gap behavior in the vor- 
tex core and no evidence for bound states. Considering 
all these points, it is fairly sure that at temperatures 
as low as 0.18K, the electronic states inside the vortex 
cores — if they exist — would not be occupied. Further- 
more, they would be localized and could not contribute 
to heat conduction. 

Extended excitations associated with supercurrents 
around vortices constitute the only other possible source 
of field-induced excess conductivity. As first pointed out 
by Volovik |juj , for an anisotropic superconducting gap, 
their contribution will dominate the variation of the den- 
sity of states, even at small magnetic fields. They are 
induced by the Doppler shift of energy spectrum due to 
a finite local value of superfluid velocity v§. This ef- 
fect is negligible in conventional superconductors at low 
magnetic fields, when the isotropic gap A >> pvs- For 
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a d x 2_ y 2 gap , on the other hand, the density of these 
excitations increases as \f~H , and simply reflects the lin- 
ear variation of the density of states with the energy, i.e. 
N(e)/Nq ~ e/A. The basic relationship can be derived 
easily : the energy of the quasi-particles being shifted 
from :£—>£ — puj(r), then the increase of the den- 
sity of states calculated on a vortex lattice cell will be 
: N(H) ~ -is §q v s {r)rdr where R ~ l/y/H. Because 

v s (r) ~ 1/r, we find that iV(if) ~ 1/R ~ Exper- 
imental evidence for the existence of this type of field- 
induced delocalized excitations in high-T c cuprates was 
first provided by specific heat measurements on YBCO 
p3| . This field-induced shift of the energy spectrum pro- 
duces effects similar to an increase in temperature. This 
equivalency is at the origin of the scaling behavior pre- 
dicted for the field-temperature dependence of specific 
heat in unconventional superconductors JlJ,|l5| . The ex- 
perimental observation of these scaling relations in spe- 
cific heat of YBCO fl6[|and thermal conductivity of UPt3 
p7[ , indicates that a convincing picture of the effects 
of the magnetic field on the density of states in uncon- 
ventional superconductors is now emerging. Recently, 
Kubert and Hirschfeld Q observed that one implication 
of this picture for the transport properties of the mixed 
state of the cuprates is to expect the thermal conductivity 
to increase with the application of a magnetic field at low 
temperatures in opposition to the decrease experimen- 
tally reported at intermediate temperatures. Our results 
constitute the experimental confirmation of such a field- 
induced increase in k by lowering temperature fis| . They 
can be naturally explained if one assumes the presence 
of nodes (or very deep minima) in the superconducting 
gap of Bi2212 at this range of temperature and field. 

In spite of a large background due to lattice contribu- 
tion to heat conduction, we can proceed to a more quan- 
titative check of this picture. The inset of Fig. || shows 
the variation of k(H)/k(0) versus \/H at T=180 mK. As 
expected, the increase in n(H) is proportional to \fH. 
However, at these low temperatures, a reliable compari- 
son with the theory would necessarily include the effects 
of the impurity band states |l9| ]. Indeed, the magnitude 
and the exact field dependence of the thermal conduc- 
tivity depends strongly on the energy dependence of the 
density of states N(e) which is in turn a function of con- 
centration of impurities and the phase shift introduced 
by their scattering potential. 

One remarkable feature of recent theoretical works on 
heat transport by field-induced quasi-particles close to 
the nodes is the prediction of a non-monotonic behavior 
of k at finite temperature [pOP] . This is a consequence of 
the energy-dependence of the relaxation time of the QP 
scattered by impurities. In presence of magnetic field, 
the Doppler shift of excitations due to the finite local su- 
perfluid velocity leads to a decrease in the relaxation time 
at low fields which can exceed the parallel rise in N(e) 
and disrupt the monotonic increase of thermal conductiv- 
ity. Hence, this theory is able to provide a unique expla- 



nation both for a monotonically increasing k(H) at low 
temperatures and a non-monotonical behavior at higher 
temperatures without invoking vortex scattering of quasi- 
particles. In this picture, the position of the minimum in 
k(H) is predicted to be H/H C2 ~ (fc^T/aAo) 2 where H C2 
is the upper critical field, Ao is the gap maximum over 
the Fermi surface, and a is a vortex-lattice-dependent 
constant of order unity. A close examination of Fig. |^ 
shows that k(T) curves at H = 600Oe and H =0 inter- 
sect at T= 0.4 K which is indicative of a temperature- 
dependent minimum in k(H). This becomes evident in 
the right panel of Fig.|[ where a logarithmic plot high- 
lights the low-field region. Assuming H C 2= 200 T, a=0.5 
and Ao = 2.14fcsT c , a rough agreement is found between 
the theoretically expected position of the minimum and 
the experimental data for T=0.55K. It is important to 
emphasize that the ordinary electron-vortex scattering 
which dominates the transport properties of the mixed 
state in conventional superconductors — which should be 
present to some extent in the cuprates — is neglected in 
this model. 
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FIG. 3. Field profile k,{H)/k,{H — 0) for various temper- 
atures. Left panel : linear- linear. Right panel : linear-log. 
At low temperatures, thermal conductivity increases with the 
magnitude of the magnetic field. H m represents the expected 
position of minimum for 0.55K(see text). The inset in the right 
panel shows k(H) / k(0) versus y/~H at 0.18 K. 

What are the implications of our results for the present 
debate on the nature of the anomaly discovered by Kr- 
ishana et ah? The grey zone in Fig. 4 represents our 
zone of exploration in the (H,T) plane where the ther- 
mal conductivity was found to increase with the mag- 
netic field due to a finite density of states. By compar- 
ing the position of this region with the extrpolation of 
the H/j (T) line to low temperatures, we can consider two 
possibilities: i) Hk(T), experimentally established from 
5 to 20 K, instead of extrapolating to low temperatures, 
increases above the grey region as shown by the dot line. 
This is improbable, ii) HuiT) extrapolates to low tem- 
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peratures with roughly the same slope(the solid line was 
theoretically predicted by Laughlin Q),but this transi- 
tion can not be seen in the thermal conductivity profile 
simply due to the lack of resolution. In this case, the 
basic augmentation of the thermal conductivity with a 
so low magnetic field imply that this new high field state 
must have a non-zero residual density of states, and con- 
sequently deep minima (or nodes) in the superconducting 
gap. 



10 



V. 


"plateau" region m. 








f/ 

LaughRrrs Line 
/ d^y state 



T (K) 



not symmetric, which would be the case in presence of 
a current associated with broken time-reversal symmetry 
in the superconducting state. Some other experimen- 
tal results have been recently interpreted as signatures 
of a splitting of the d x 2_ y 2 state. For example, Franz 
and Tesanovic |j| have proposed that the existence of a 
d x 2_ y 2 + id xy state could be at the origin of the bound 
states observed in the vortex cores of YBCO by tunnel- 
ing spectroscopy JT0| ] . Furthermore, the observation of 
a spontaneous split of the zero-bias conductivity peak at 
7K in YBCO pi] , was attributed to a phase transition at 
the sample surface from a pure d-wave to a time-reversal 
violating state. 

In conclusion, our results show that the new high-field 
superconducting state in Bi2212 identified by the emer- 
gence of a plateau in the field-dependence of the thermal 
conductivity cannot be simply attributed to a vanishing 
of low-energy excitations in the superconducting state. 
At low-temperature, k increases with a weak magnetic 
field, implying a non-zero density of states. 

We are indebted to L. Fruchter for magnetization mea- 
surements. Fruitful discussions with him and with L. 
Taillefer, P. Hirschfeld and Y. Barash are acknowledged. 
This work was partly supported by Grant-in- Aid for Sci- 
entific Research from the Ministry of Education, Science, 
Sports and Culture, Japan. 



FIG. 4. k was found to increase with magnetic field in the 
grey zone of the (H,T) plane. Full circles represent experi- 
mental Hfc(T) extracted from Ref. 1. The solid line is the 
theoretical prediction by Laughlin. 



These observations do not rule out a phase transition. 
Indeed, the kink observed is robust from one sample to 
another, and is not obviously related to a phase transi- 
tion in the vortex lattice. However, if this phase transi- 
tion in the electronic system exists, its nature or conse- 
quences on the electronic heat carriers are different from 
what is simply expected for a state with a large uniform 
gap suppressing the excitations of the superconducting 
state in the whole sample. Already, when presenting his 
model on field- induced transition to a d x 2_ y 2 + id xy || 
state, Laughlin pointed out that the expected gap is too 
small to account for a complete vanishing of QP trans- 
port. Moreover, he stressed the inherent difficulty in ex- 
plaining a large field-induced gap exceeding ksT c , the 
temperature at which the transition occurs. 

At this stage, a broken time-reversal symmetry associ- 
ated with the splitting of d x 2_ y 2 state remains the most 
serious track to explain the origin of the strong hysteretic 
behavior. Indeed, the field profile in the sample is asso- 
ciated with a superfluid current of the order of J c , the 
critical current. When reversing the magnetic field ramp- 
ing direction, the most dramatic change in the sample is 
the direction of the field profile and the way of circula- 
tion of the screenings currents. The sensitivity of ther- 
mal conductivity to the direction of circulation of these 
screenings currents may indicate that these situations are 
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